Dragging a polymer chain into a nanotube and subsequent release 
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We present a scaling theory and Monte Carlo (MC) simulation results for a flexible polymer chain 
slowly dragged by one end into a nanotube. We also describe the situation when the completely 
confined chain is released and gradually leaves the tube. MC simulations were performed for a 
self-avoiding lattice model with a biased chain growth algorithm, the pruned- enriched Rosenbluth 
method (PERM). The nanotube is a long channel opened at one end and its diameter D is much 
smaller than the size of the polymer coil in solution. We analyze the following characteristics as 
functions of the chain end position x inside the tube: the free energy of confinement, the average 
end-to-end distance, the average number of segments imprisoned in the tube, and the average 
stretching of the confined part of the chain for various values of D and for the number of repeat 
units in the chain, N. We show that when the chain end is dragged by a certain critical distance 
x* into the tube, the polymer undergoes a first-order phase transition whereby the remaining free 
tail is abruptly sucked into the tube. This is accompanied by jumps in the average size, the number 
of imprisoned segments, and in the average stretching parameter. The critical distance scales as 
x* ~ ND 1-1 '". The transition takes place when approximately 3/4 of the chain units are dragged 
into the tube. The theory presented is based on constructing the Landau free energy as a function 
of an order parameter that provides a complete description of equilibrium and metastable states. 
We argue that if the trapped chain is released with all monomers allowed to fluctuate, the reverse 
process in which the chain leaves the confinement occurs smoothly without any jumps. Finally, we 
apply the theory to estimate the lifetime of confined DNA in metastable states in nanotubes. 



I. INTRODUCTION 

Recent developments in fabrication of nanoscale de- 
vices and in single-chain manipulation techniques open 
possibilities for a broad range of applications in biotech- 
nology and materials science [H,[3,[|l,i,[i,[3. In partic- 
ular, well calibrated nanochannels were produced in fused 
silica substrates by lithography methods with the widths 
in the range of 30 to 400 nm, which were used to study 
the confinement of single A-phagc DNA molecules driven 
electrophoretically into these nanochannels 0] . The per- 
sistence length of DNA under conditions used in these 
experiments is about 50 nm while its contour length was 
about 1000 times larger. This means that except for the 
case of the narrowest channels DNA behaved essentially 
as a long flexible macromolecule on the relevant length 
scales. The aim of this paper is to elucidate the subtle 
physics behind the process when a long flexible chain is 
slowly dragged by one end into a nanochannel. We also 
show that this process is qualitatively different from what 
happens when a confined chain is released and leaves the 
nanochannel by spontaneous thermal motion. 

Our approach is based on using the most general re- 
sults of the scaling theory that neglect the small-scale 
details of the system under consideration, namely the 
particular form of the interaction potentials, the flexibil- 
ity mechanisms, etc. It is well known that the scaling ap- 



proach does not allow to calculate non-universal numer- 
ical coefficients that are model-dependent. To verify the 
prediction of the analytical theory we have carried out de- 
tailed Monte Carlo (MC) simulation by u sing th e pruned- 
enriched Rosenbluth method (PERM) HElME]- Thc 
simulations are based on the simplest model of polymers, 
namely self-avoiding walks on a cubic lattice. Fortu- 
nately, the general ideas of scaling guarantee a universal 
behavior within broad limits of parameters. 

The properties of a single macromolecule confined 
in a tube have been studied extensively for decades, 
both by analytical theory and by numerical simula- 
tions for various models of flexible and semi-flexible 
chains [HI, [3 [ID, fig . For a homogeneous confined 
state there are scaling predictions [T^| concerning various 
chain characteristics which were tested by MC simula- 
tions. Our main interest here is in the non-homogeneous 
flower-like states where the confined part of the chain 
inside the tube forms a stretched stem and the free tail 
still in solution forms a coiled crown. This type of con- 
formations appears in a variety of situations including 
translocation through a thick membrane [18| as well as 
the escape transition produced by c omp ressin g a grafted 
chain between flat pistons [ijj [2(1 l2ll I22I [231 ] . One of 
the goals of this paper is to demonstrate that the role 
of these conformations in the situations mentioned above 
are quite different depending on whether one of the chain 
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FIG. 1: Schematic drawings of a flexible polymer chain with one end dragged into a nanotube in a quasi-equilibrium process. 
The fixed chain end characterized by its coordinate x experiences an reaction force / that balances the pull of the undeformed 
swollen coil outside the tube (a)(b). At the transition point x* , the remaining tail is sucked into the tube abruptly by a uniform 
shrinking of all blobs in the stem (b), (c), and the reaction force becomes zero, / = 0. At the transition point conformations 
(b) and (c) coexist. As long as the chain is fully confined in the tube, no reaction force appears at the fixed chain end (d). 



ends is fixed in the confinement region or not. 

The paper is organized as follows: We start by pre- 
senting the main qualitative results in a series of simple 
pictures visualizing the conformational changes in the 
process of dragging the chain end slowly into the tube 
and its de-confinement upon subsequent release. In Sec- 
tion III we give the results of the MC simulations for 
completely confined states and compare them with the 
well-known scaling prediction. Section IV describes the 
simulation results characterizing the phase transition in- 
duced by changing the chain end position inside the tube. 
In Section V we provide a theoretical description of the 
phenomenon based on the Landau free energy approach 
and compare it to the MC results. The process of spon- 
taneous de-confincment is analyzed in Section VI and 
followed by a general discussion. 



II. QUALITATIVE PICTURE 

Figure [T] presents a sequence of chain conformations 
when the position of one of the chain ends is progres- 
sively moved quasi-statically further inside the nanotube. 
At each moment, this fixed chain end is characterized by 
its coordinate x which is controlled externally. The tail 
outside the tube is a practically undeformed swollen coil 
while the part of the chain inside forms a one-dimensional 



string of blobs. The chain can be presented as a flower 
structure with one-dimensional structure of the stem. 
The number of confined (imprisoned) chain units grows 
linearly with x, at least in the initial stages of the pro- 
cess. The fixed chain end experiences a force which is 
due to the tail still not being confined. To counterbal- 
ance this force, a reaction force / directed into the tube 
will appear. Once the chain is confined completely (Fig- 
ures [Tb and [TJi) there is no tail outside the tube and 
the reaction force at the controlled chain end disappears. 
In this state, the chain is homogeneously stretched and 
its stretching degree is due solely to the confinement ef- 
fect pjj ■ It is clear that the stem of the partially confined 
conformation is stretched more strongly since the confine- 
ment effects arc augmented by the additional stretching 
force (this effect is symbolically indicated by the defor- 
mation in blob shape shown in the Figures QJi andO})). 
The difference in the deformation free energy leads to a 
phase transition - an abrupt "slurping" of the remaining 
tail accompanied by a uniform shrinking of all blobs in 
the stem. (Of course, in a strict sense true phase transi- 
tions can occur only in the thermodynamic limit, which 
would require that both the length of the chain and the 
length of the pore are infinitely large; however, as we 
shall see, the rounding of the phase transition caused by 
finite chain length is not too significant.) As a result the 
length along the tube of the completely confined chain is 
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FIG. 2: Schematic drawings of a released flexible polymer chain. Since all the chain units are free to fluctuate, no reaction 
force is needed to counterbalance the force due to the tail outside the tube. The stretching of the stem is only dependent on 
the confinement effect at any moment, (a) and (b). 



less than the length of the stem just before the transi- 
tion, as illustrated in Figures [TJd and [TJ:. Note that the 
process described above has to be distinguished from the 
case when the chain is driven into the pore by applying 
a constant force, as in electrophoresis. 

The next figure (Figure [2]) illustrates the de- 
confincmcnt of a released chain. All the chain units are 
free to fluctuate and there is no reaction force. There- 
fore, the stem stretching is only due to the confinement 
effect at any particular moment. The gradual decrease 
in the number of imprisoned units is not accompanied by 
any jumps. 



III. FULLY CONFINED CHAIN: 
AND SCALING 



MC RESULTS 



The scaling picture of a fully confined chain inside a 
tube of diameter D is very simple. It is formed by a 
string of non-overlapping blobs of size ~ D, each blob 
containing g monomer units [171 ] . The size of the blob is 
related to g by 



D~ag" 



(1) 



where a is the length of a monomer unit, and v = 
0.58765(20) is the scaling index for a three-dimensional 
self-avoiding chain [ll|. Let's define the number of blobs 
by relation 



n b = N/g = N{D/ a y 1/lJ 



(2) 



In this case the average end-to-end distance of the chain 
in a fully confined (imprisoned) state is 



(3) 



where A lrap is a model-dependent numerical coefficient. 
In eqs [T][3] we have assumed that the pore diameter D 



is large enough, so that the number of monomers inside 
a single blob is large enough so that the scaling relation 
cqQ] holds; further more correction terms to the scaling 
relations are omitted throughout. 

A stretching parameter S- lmp — Ri mp /Na describes the 
average stretching of the chain in an imprisoned state. 
Since the end-to-end distance of the fully imprisoned 
chain is proportional to A, the stretching parameter is a 
function of the tube diameter only and is given by 
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A imv {D/a) l - l /» 



(4) 



For a fully imprisoned chain the free energy of confine- 
ment per blob scales as ksT. Thus, using eq[2]the free 
energy of a fully imprisoned chain is 



(5) 



where i?i mp is another numerical coefficient 



The factor 

of ksT is absorbed in the free energy throughout the 
paper hereafter. 

For our simulations, chain lengths are up to 44000, 
tube diameters are up to D = 97, and a = 1 which is the 
lattice spacing. Simulation data for the rescaled average 
root mean square (rms) end-to-end distance Rim P /(ni,D) 
are displayed in Figure [3] depending on the blob number 
rib- It is evident from the Figure [3] that the normal- 
ized chain size reaches a practically constant value when 
the chain contains more than two blobs. The limiting 
constant value of Ri mp /(ni,D) is just the numerical co- 
efficient Ai mp = 0.92 ± 0.03 in eq [3] As noted above, 
the scaling description is expected to become exact in 
the asymptotic limit D — > oo, and hence in the consid- 
ered range of not very large D a small systematic depen- 
dence of the curves in Figure [3^ is evident, leading to 
the spread of values of the coefficient A- imp , as expressed 
by the quoted uncertainty. Figure [3] shows that at values 
nt, < 1, the unperturbed coil size smaller than the tube 
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FIG. 3: (a) Rescaled average end-to-end distance Ri mp /(nbD) 
against the number of blobs rib = N(D/a)~ 1 ^ L ' for fully con- 
fined chains. The thick solid line with arrow is J?i mp /(nf,D) = 
Aimp with Amp = 0.92(3). (b) The log- log plot of the same 
data as (a) shows that for rib > 2 the data collapse and allows 
to estimate j4i mp very accurately. The dotted straight line 
with slope v — 1 shows the scaling law in the regime of wide 
tubes with D > i?i mp . 



diameter, there is another scaling regime of weak confine- 
ment where i?i mp ~ N v . Although the crossover between 
these two regimes is of general theoretical interest, we are 
not concerned with it in this paper. 

The rescaled free energy F; mp (counted from the ref- 
erence state of a self-avoiding coil in free space) for the 
fully imprisoned chain is presented in Figure 2] vs. the 
number of blobs rih. The limiting constant value gives 
the coefficient B imp — 5.33 ± 0.08 in eq which rep- 
resents the free energy per blob in UbT units. Once 
again we see the other scaling regime of weak confine- 
ment if the number of blobs in the chain is less then one 
blob. Summarizing the presented data we can conclude 
that using the definition of the number of blobs given 
by eqd]the size of the blob in our model is close to D, 
the end-to-end distance Ri mp ~ 0.92nf,D, and the free 
energy per blob F imp /nb ~ 5.33 is close to 5 fceT. Com- 
paring our result for the free energy per blob, -Ftubc/^b 
(-Ftube = -fimp) with that of chains confined in a slit on a 
lattice [ll|, Fgut/rib ~ 2.10, and a off-lattice model [26| . 



(b) 



E 
LL 




0.1 1 

n b = N(D/a)~ 



1/v 



FIG. 4: (a) Rescaled average free energy Fimp/nj, for fully 
confined chains against the number of blobs rib ■ The thick 
solid line with arrow is Fimp/rib = -Bimp with Bi mp = 5.33(5). 
(b) Same as (a), but on a log- log plot to show that for n;, > 1 
the data allows to estimate -Bi mp very accurately. 



Fgiit/nb ~ 2.03, we find that F tu b c ~ 2-F s iit. It is not 
surprising because a polymer chain confined in a slit is 
compressed in one direction and in a tube in two direc- 
tions. 



IV. PHASE TRANSITION: EQUILIBRIUM 
CHARACTERISTICS 

For our simulations, single polymer chains are dragged 
into a tube with diameters D — 17, 21, 25, and 29, and 
the chain length is up to TV = 17000. Simulation re- 
sults of the free energy relative to a self-avoiding coil are 
presented in Figure [5] as a function of the end monomer 
position inside the tube, x, normalized by the tube diam- 
eter D. It is clear that there are two branches of the free 
energy. Initially, the free energy increases linearly with 
x as more and more blobs are driven into the tube. The 
slope of the free energy as a function of x has the meaning 
of the average reaction force acting on the end monomer. 
Deviations from the linearity near the origin occur when 
only the number of blobs inside the tube is of order one 
or less. At large enough values of x all monomeric units 
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FIG. 5: Free energy of chains with chain lengths JV = 500 and 
N = 1500 relative to a self-avoiding coil, F(N,x,D), plotted 
against x/D for tube diameters D = 17, 21 and 25. The solid 
straight line is F$ = 4.23a;/U and gives the best fit of the data 
of chains in a flower state. All horizontal lines indicate the 
values of Fi mp for chains in an imprisoned state at fixed chain 
length N. The intersections of the solid line and horizontal 
lines indicate the transition points x* / D. 
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FIG. 6: Transition points x* as obtained from Figures(5] plot- 
ted against N(D/a) 1 ~ 1 ^ l/ for various values of N and D. The 
dotted line has slope Bi mp /Bfl = 1.25 . 



are confined and the free energy does not depend on the 
position x any more. An abrupt change in the slope of 
the free energy indicates a first-order transition. The lin- 
ear branch describes a partially confined "flower" state. 
On this branch the data points for different values of N 
and D collapse onto the same universal curve in the cho- 
sen coordinates. The x-independent branch corresponds 
to a completely confined state discussed above in Section 
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FIG. 7: (a) Average fraction of imprisoned units, Ni mp /N, 
plotted against the reduced end coordinate x/x* for rib ~ 60 
and for various values of D. The solid line gives the theoretical 
prediction that the fraction increases linearly with x below the 
transition point, and it is in perfect agreement with our data. 
At the transition point x/x* = 1, Ni mp /N jumps up from 0.76 
to 1, and the relative reduction in the number of imprisoned 
monomers is then Ajv ~ 0.24. (b) JVi mp /JV vs. x/x* the 
transition region, x/x* — 1, for D = 21 and for different 
number of rib, displaying the rounding of the transition is due 
to the finiteness of the number of blobs. 



III. Altogether the results can be summarized as follows: 



F(N,x,D) = 



_ J B iTav N{D/a)-V v =F h 



imprisoned state 
Br(x/D) — Fa flower state 

(6) 

where £?i mp = 5.33(8) and Bq = 4.23(6) were obtained 
from Figure 2] and Figure [5] Physically, both formulas 
state that the free energy is proportional to the number 
of blobs inside the tube. However, we would like to point 
out that the numerical coefficients are different, and the 
most important source of this difference is due to the 
extra stretching of the flower stem as compared to the 
relaxed fully confined state. 

The intersection of the two free energy branches defines 
the transition point, or the critical distance x* away from 
the open end of the tube. Equating the two expressions 
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FIG. 8: (a) Reduced rms end-to-end distance R/Rim P , plot- 
ted against the reduced end coordinate x/x* for rib = 60 and 
for various values of D. Ri mp is the average rms end-to-end 
distance for a fully confined chain. Near x = 0, chains be- 
have as a self-avoiding coil, R/Rim P depends on both N and 
D. Approaching the transition point from below, a linear be- 
havior appears according to the scaling prediction. At the 
transition point x/x* = 1, the value of R/Ri mp jumps down 
from 1.38 to 1 and the reduced jump of end-to-end distance 
A R = (1.38 - 1)/1.38 w 0.27. (b) R/R imp vs. x/x* around 
the transition region, x/x* = 1, for D = 21 and for different 
number of m. 
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FIG. 9: (a) Reduced average stretching parameter < s > 
/Simp, plotted against the reduced end coordinate x/x* for 
rib = 60 and for various values of D. S mp is the average 
stretching of a fully confined chain, and it does not depend 
on x/x* . The average stretching of the chain in a flower state 
is the stretching of the stem and it is also independent of x/x* 
except the region that near x = where it depends on both 
N and D. At the transition point x/x* — 1, the value < s > 
/Si mp jumps down from 1.75 to 1 and the reduced jump of 
the stretching parameter is As = (1.75 — 1)/1.75 ~ 0.43. (b) 
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x/x* around the transition region, x/x* = 1, 



for D — 21 and for different number of blobs rib. 



of eq[6]one obtains: 

x* = (S imp /5 fl )iV(D/a) 1 - 1 /- = 1.26(4)iV(/J/ a ) 1 - 1 /- 

(7) 

MC data for the transition points displayed in Figure [B] 
are in full agreement with eq[7] 

The average fraction of imprisoned units, N^p/N, as 
a function of the reduced end coordinate x/x* is shown 
in Figure [7J It is clear that this fraction increases lin- 
early with x and jumps up to 1 at the transition point. 
This jump represents the abrupt "slurping" of the tail 
which constitutes approximately 1/4 of the total num- 
ber of units, N . As mentioned above, truly sharp jumps 
can occur in the thermodynamic limit [N — * oo) only, so 
when one looks at the behavior of Ni mp /N with high res- 
olution on the scale x/x* for various sizes of the system, 
one can resolve the finite-size rounding (Figure Uh)- The 



next graph (Figure [8]) represents the change in the rms 
end-to-end distance reduced by the corresponding value 
for the completely confined chain, R/ R lmv . The starting 
value at x — represents an unconfined coil with R w N u 
so that the reduced value depends on both N and D, i.e., 
R/Rimp oc N u /{ND 1 - 1 /"). Then, a linear dependence 
appears which is terminated at the transition point. The 
average size jumps down by approximately 25%. Finally, 
the average stretching parameter < s > of the confined 
part of the chain is shown in Figure [5] as a function of 
the same reduced distance x/x*. (For the precise defini- 
tion of the stretching parameter s, see the next section). 
It is reduced by the value S lmp characterizing the fully 
confined chain, independent of x and given by eq 2] In 
a flower state the stretching parameter Sr describes the 
average stretching of the stem. It is independent of x 



7 



also and is about 1.75 times larger than that for the fully 
confined chain. A small deviation exists near x = only. 
Beyond the transition point Sr jumps down to the 5i mp 
value. 



V. LANDAU THEORY 

Unlike the case of critical phenomena, the Landau the- 
ory approach is very well suited for analyzing first-order 
transitions, including possible metastable states. The 
idea is to first subdivide all configurations into subsets as- 
sociated with a given value of an appropriately chosen or- 
der parameter s that allows to distinguish between differ- 
ent states or phases. Landau free energy $(s) is the free 
energy of a given subset, and is therefore a function of the 
order parameter. The minimum value of the Landau free 
energy is attained for the subset that contains most of the 
equilibrium configurations, and therefore coincides with 
the equilibrium free energy of the system, F. Far enough 
from the first order transition point, the Landau free en- 
ergy has only one minimum. However, near the transition 
the function is expected to have two minima (the deeper 
one is stable and the other is metastable). Exactly at 
the transition point, both minima are of equal depth. 
We define the order parameter as the chain stretching 
in the fully confined state, s = r/(N- lmp a) (iVi mp = N) 
where r is the instantaneous end-to-end distance of the 
chain, or as the stretching of the stem only in the flower 
state, s = x I (rid) where n is the instantaneous number 
of confined monomers in the stem, and x is the length of 
the stem. The Landau function consists of two branches 
that have to be introduced separately. For fully confined 
configurations the Landau free energy, up to an additive 
constant, is directly expressed in terms of the distribu- 
tion of the end-to-end distance P(r\N, D) of a chain in 
the tube: 

$ imp {r/Na) = const - \nP(r\N,D) . (8) 

There exists no closed formula for such a distribution of 
confined chains with excluded volume interactions. In 
Ref. [HI the (non-normalized) distribution for the gy- 
ration radius was studied numerically and the following 
scaling form was proposed: 

]nP(R g \N,D) = -N{D/a)-^ v A[ u - a + Bu s ] (9) 

where a = (3v — 1) _1 ,6 = (1 — v)^ 1 and u — 
{R g /Na)(D/a)- 1+1 / v . The parameters A and B are non- 
universal numbers of order unity and they do not depend 
on N or D. The first term describes the concentration 
effects in the des Cloizeaux [24| form, and the second 
term is the Pincus [25[ scaling form of the stretching free 
energy. From our MC simulations, we have found that 
the end-to-end distribution P(r\N,D) and the equilib- 
rium free energy F are well described by a similar scaling 
formula corrected by an additional r-independent term, 



namely the Landau free energy for the imprisoned state 
is then 

$ imp (s) = N{D/a)- l ' v A [u- a + Bu s + C] (10) 

where u is now related to our order parameter s 

u = {r/Na){D/a)- 1+1 / v = s{D/a)- 1+1 ' v . (11) 

This branch is limited to the range of values for the or- 
der parameter, < s < x/Na. In the thermodynamic 
limit, the average value of the order parameter of a fully 
confined chain, Si mp =< s >, is found by locating the 
minimum of $i mp (s), i.e., d$ lmp (s) / ds — at s = S- lmp , 
and the corresponding minimum of the Landau free en- 
ergy is the equilibrium free energy F. Using eq 1 1 H we 
obtain the equilibrium average value of the end-to-end 
distance i?i mp —< r >, 

Rim P = N <s>=N{D/a) 1 - 1 / v u 1 (12) 

where u\ = {a,/8B) l /( a+ ^ gives the position of the mini- 
mum of the function f(u) — u~ a + Bu 6 + C . Comparing 
this result with the MC data of D = 17, i.e. eq[3]with 
taking Ai mp = 0.94, we immediately determine the nu- 
merical value of B = 0.67. Using eg 1101 the equilibrium 
free energy Fi mp is the Landau free energy at u = U\: 

fimp = N(D/a)- 1 ' v A (1.67 + C) (13) 

Comparing this again with the MC data of D — 17, i.e. 
eq[5]with taking i?i mp = 5.38, we get a relationship be- 
tween the coefficients A and C . The last condition that 
eventually fixes all the numerical coefficients A, B and 
C of the Landau function is obtained by analyzing its 
second branch. 

For the partially confined chains in the flower state, 
since in fact only n monomers that comprise the stem of 
the flower give contributions to the free energy, instead 
of N and r, we use n and x in eq 1101 The formula of the 
Landau free energy is therefore, 

$ fl (s) = ^A [u-*- 1 + Bu 5 - 1 + Cu- 1 } , fl>^(14) 

where u is given by 

u = {x/na){D/a)- 1+1/l/ = s(D/a)- 1+l ' v . (15) 

The average value of the order parameter Sq, and the 
equilibrium free energy Fr for the flower state are ob- 
tained by the same procedure as that for the imprisoned 
state. Demanding that Sg and the coefficient with the 
x/D factor in eq[T3]both coincide with the MC results, 
we are uniquely fixing the numerical values of parameters 
A = 1.48, B = 0.67 and C = 1.98. 

VI. COMPARISON WITH MC RESULTS. 

For our simulations, the Landau free energy as a func- 
tion of s, $(iV, x, D, s), is given by taking the logarithm 
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FIG. 10: Landau free energy divided by N, $(N,x, D, s)/N, 
plotted against the order parameter s for various values of 
x/x* and for D = 17 (a) and D = 21 (b). The predicted 
Landau free energy function 3>i m p, eg 1101 in the imprisoned 
regime, and $a, eq 1141 in the flower regime, with A — 1.48, 
B = 0.67, and C = 0.98, give a good fit to the MC data. 



of the properly normalized accumulated histogram of 
stretching parameter s. Results for tubes with diam- 
eters D = 17 and D — 21 near the transition point 
x*/Na ~ l^D/a) 1 - 1 /" {eq[7]} are shown in FigureOU 
Two branches of the analytical Landau function given 
by eqs [TO]and[EI] for A = 1.48, B = 0.67, and C = 1.98 
are also presented in Figure [TOJ We see that in both 
the flower and the imprisoned states, MC data are in 
perfect agreement with the theoretical predictions for a 
wide range of s at any x, and D although the coefficients 
A, B, and C are determined by using only the MC es- 
timates of average stretching parameters S- imp and Sa , 
and of average free energies F; mp and Fa for D — 17. 

From the analytical Landau function in the imprisoned 
state and in the flower state, eqs[TUland[T4"l and the de- 
termined values of the coefficients A, B, and C, we obtain 
the following scaling relationships of the equilibrium or- 
der parameters s 



OMiD/a) 1 - 1 ^ = S imp 
lMiD/a) 1 - 1 ^ = Sb 



imprisoned state 
flower state 

(16) 



and the equilibrium free energies F, 

_ ( 5A0N(D/a)-^ v = F imp imprisoned state 
eq ~ \ A.27x/D = F & flower state 

(17) 

The transition point is found from the condition that the 
two minima of the Landau free energy function are of 
equal depth. Using eqs[T7|we get 



- 1.26(D/a) 1 - 1 ^ 
Na v 1 ' 



(18) 



The reduced jump of the order parameter is therefore 



Ac 



0.43 



(19) 



The average number of units dragged into the tube for 
an imprisoned state is A^ mp = N, while for the coex- 
isting flower state we have only A^ mp =< n >= x/Sr 
monomers. Using eqs [Tf)] and [TS] we obtain the relative 
reduction in the number of imprisoned monomers 



N - x/S a 
N 



0.23 



(20) 



Finally, the reduced jump of the end-to-end distance is 
obtained by using eqs [16] and [T8l 



R = 



0.25 



(21) 



Equations [1911211 show that the sizes of jumps in S, 
-Wimp/-W> an d R are universal quantities (i.e., indepen- 
dent of D and N), and they satisfy the following relation 



(1 - A R )(1 - A N ) = (1 - A s ) 



(22) 



Comparing with the MC results shown in Figures 091 we 
see that the Landau theory gives a good qualitative and 
quantitative agreement. 



VII. METASTABLE REGIONS AND SPINODAL 
POINTS 

When the chain is dragged into the tube by one end 
slowly (quasi-statically), the number of imprisoned units, 
Mmp as a function of x grows linearly up to the transi- 
tion point x* and then jumps to N corresponding to a 
full confinement as shown in Figure [7J This would mean 
that at any value of x a complete equilibrium is achieved 
and only the lowest minimum of the Landau free energy 
is populated. However, Figure [TTk shows clearly that for 
x > x* there still exists another minimum of the Landau 
free energy representing the metastable flower state. The 
barrier height per blob turns out to be a universal func- 
tion of the reduced coordinate x/x* that is displayed in 
Figure [Tib. So, if the chain end is moved into the tube 
relatively quickly compared to the metastable lifetime, 
the chain will be trapped in the flower state, and the 
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0.7 0.8 0.9 1 1.1 1.2 1.3 

Xsp/X x/x* Xsp/X* 

FIG. 11: (a) Landau free energy per segment, 
<&(N,x,D, s)/N, versus order parameter s close to the 
transition point for x/x* = 1.05 and for D — 17. The 
two branches intersect at some intermediate state with 
s = Sinter- From the difference between &(N,x,D,s)/N 
at the intersection point and at the higher minimum one 
finds U/N = 0.0021. (b)Barrier heights per blob, U/n b , 
estimated as described in (a), vs. x/x*. In the coordinates 
used the curve is universal (independent on D, and TV) At 
the transition point x/x* = 1 the barrier height is maximal. 
At the two spinodal points x$ /x* = 0.78 and xip 1 /x* = 1.33 
the barrier vanishes. 



X 



0.1 0.2 0.3 0.4 0.5 

s 

FIG. 12: (a) Number of imprisoned units (normalized to 
unity), Ni m p/N, plotted against the coordinate of the end 
monomer position of a dragged chain, x. If the chain end is 
moved into the tube very slowly, Nim P /N of chains in equilib- 
rium states follows the solid curve. If the chain end is moved 
quickly, the chain may be trapped in the metastable states 
(dashed curves) and a hysteresis loop appears, (b) Landau 
free energy divided by N, <&(N,x,D, s)/N, as a function of 
the order parameter s for D = 17 at the transition point x* 
and at the two spinodal points xip and xip . In order to dis- 
tinguish these three curves well, a constant value A is added 
to $(N,x,D,s)/N. A = 0.02, 0, and -0.02 for x = x$ , x* , 

(2) 

and x = xip , respectively. 



number of imprisoned units will keep increasing linearly 
as shown by the dashed line in Figure [L2r . As the tail 
decreases, so does the barrier height, as demonstrated in 
Figure [Tib until the metastability is completely lost at a 

(2) . (2) 

spinodal point x = x sp . The spinodal value Xs P is de- 
fined by the condition that the flower minimum coincides 
with the matching point (Figure [12b) leading to 

zg) = NaSfi (23) 

where <Sfl is given by eq 1161 

If the chain is in a fully confined state and its end 
is moved back to the tube entrance quasi-statically, the 
equilibrium curve is retraced. On the other hand, a 
metastable imprisoned state also appears at x < x* . Its 
local stability is lost at the other spinodal point, 

z« = NaS imp (24) 



where Si mp is again given by eqUBJ It follows from the 
physical meaning of S- unp that xip coincides with the 
equilibrium end-to-end distance for a fully confined chain. 
Note that both spinodal points scale in the same way 
with N and D, x sp - ND 1 ^ 1 / U although with different 
numerical prefactors. The hysteresis loop associated with 
the metastable states is displayed in Figure [T2"a. Let us 
estimate the lifetime of a metastable state of a single A 
phage DNA confined in a nanochannel 0] with the follow- 
ing parameters: contour length L = 16/im, persistence 
length a = 50nm, tube diameter D — lbOnm. This gives 
the number of blobs {L/ a){D / a)- 1 /" = 50. The lifetime 
of a metastable state (mean first passage time) can be es- 
timated as r ms = tq exp(L T /rCBT) where U is the height 
of the barrier separating the metastable minima and the 
interaction point. A characteristic relaxation time tq for 
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the DNA molecule estimated from autocorrelated exten- 
sion fluctuations is close to 1 sec. The barrier height 
near the transition point is about 0.38 ksT per blob ac- 
cording to Figure [TTb . With these parameters an esti- 
mate for the lifetime of a mctastablc state is astronomical 
i T ms ~ 10 13 s, leaving no chance to observe the equi- 
librium transition experimentally and making hysteresis 
effects inevitable. However, since the number of blobs 
depends strongly on the width of the tube one can ex- 
pect that experiments with a wider tube would be much 
closer to equilibrium. As an example, for the same DNA 
molecule in a tube with D = 500nm a similar estimate 
gives nfc = 10 and r ms ~ lmin. 



VIII. ESCAPE OF A RELEASED CHAIN 



M.™ ensemble 




x ensemble 



<x> 



FIG. 13: The number of imprisoned units (normalized to 
unity), Nimp/N, plotted against the end monomer position 
of the released chain, x, in two ensembles. In x-ensemble, the 
end position changes quasi-statically. In A^mp-ensemble, the 
number of imprisoned units changes quasi-statically. 



We are now in a position to address the second part 
of the problem announced in the title of the paper. We 
have demonstrated that dragging a chain into a tube by 
its end involves a first order transition accompanied by a 
jump-wise change in the chain conformation. The ques- 
tion is then, whether its escape back from the tube upon 
release will also involve a jumpwise transition. To clar- 
ify the situation we recall that we were trying to simu- 
late and describe theoretically an experimental situation 
where the position of one chain end serves as a parameter 
controlled by external means, e.g. by using optical tweez- 
ers. Theoretically, this implies a statistical description in 
a constant x ensemble in which statistical averaging is 
done over all internal degrees of freedom at fixed val- 
ues of x (together with other parameters such as N, D, 
and temperature T). The averaged quantities, e.g. the 
average number of imprisoned units, are thus functions 
of x. Their variations with x describe the response of 
thermally fluctuating variables to a change in external 
parameters. 

In a gradual escape of a released chain initially confined 
in a tube the coordinates of all segments, including both 
ends, are themselves subject to thermal fluctuations. An 
important question to be addressed is: what is the appro- 
priate statistical ensemble to describe this process? Let 
us first assume that the position of the distant end inside 
the tube, x , is indeed a dynamic variable that changes 
much slower than the other degrees of freedom. (We 
argue below that this assumption is generally incorrect 
unless specific mechanisms are in place to achieve this 
effect). Then the evolution of x itself will be governed by 
the free energy profile F(x) discussed earlier and shown 
in Figure El Dynamically, the chain will diffuse along 
the flat horizontal portion of the slope and slide down 
the slope towards the de-confined state at x < 0. In the 
process, both the average x, < x >, and the average n, 
Ni m p, will be changing with time but iVi mp as a function 
of < x > defined parametrically will follow the retraced 
equilibrium curve in Figure [T5J 

Now the relaxation time for the number of imprisoned 
units, iVj mp , presumed to be governed by faster dynam- 



ics is to be estimated. For x in the range between two 
spinodal points xi]} < x < xip^ , there are two minima 
separated by a barrier, see Figure W2h leading to a very 
slow relaxation that involves barrier crossing. It is clear 
that this contradicts the assumption that N- lmp adjusts 
quickly to any change in x unless the motion in the x co- 
ordinate is specifically slowed down by some additional 
mechanism. As an example of such a mechanism, one 
could envisage a situation where the distant chain end is 
modified to have a sticky anchor attached. This would 
not affect the fluctuations of the end nearest to the tube 
opening and thus would not slow down the process of 
equilibrating the number of imprisoned segments that 
involves expulsion of a tail. A relatively large colloidal 
particle attached to the distant chain end for the purpose 
of using optical tweezers may produce a similar effect. 

A more appropriate ensemble seems to be the one 
where the number of imprisoned segments is assumed 
to change quasi-statically while the end position x is ad- 
justed by thermal fluctuations. In the absence of any 
external force the relaxation of x is purely diffusive if 
the chain far inside the tube with A^ mp = N. Once the 
process of gradual escape starts and N\ mp < N the co- 
ordinate x relaxes to produce the equilibrium stretching 
of the remaining confined part at the value s = S{ mp . 
This relaxation is never controlled by barrier crossing 
rate. The evolution of < Ni mp > (t) itself can be pic- 
tured as a slide down the linear slope of the free energy 
F{Ni mp ) — Ni mp Bi mp (D I 'a)~ 1 ' v towards the minimum 
n = 0. The parametric dependence of < N{ mp > vs. 
< x > in this process is also shown in Figure [T3l 

The use of the fixed i\Zi mp ensemble would be clearly 
justified if the evolution of the variable N- imp was to be ex- 
plicitly slowed down without affecting the relaxation rate 
of the distant end. An experimental situation where this 
slowing down could be realized if the chain was escaping 
from the tube trough a partially blocked opening, similar 
to a setting normally assumed in a problem of transloca- 
tion through a narrow hole in a thick membrane. 



11 



It follows from the above discussion that the behavior 
of a released chain escaping from a tube is a problem 
of real polymer dynamics which is far from being well 
understood. By using quasi-equilibrium statistical en- 
sembles we were able to clarify the two limiting dynamic 
cases when one of the global variables is much slower 
than the other. Although both limiting results may be 
applicable to real experimental situation provided some 
additional modifications of the basic setting are intro- 
duced, one could speculate that the process of an escape 
from a tube in the simplest setting is somewhere in be- 
tween these limits. A more general and powerful ap- 
proach where both N- lmp and x are treated dynamically 
is sketched in the Appendix. 



end at a particular position x is removed, so that no 
external force ever acts at the chain end inside the tube: 
the chain then can reduce its free energy continuously by 
"escaping out" of the tube. The dynamics of this chain 
expulsion process from a tube is an interesting problem 
for further study, however. 

Another interesting extension would concern the be- 
havior of a chain dragged into a tube with attractive 
walls. We expect that such a situation could be of in- 
terest in the context of polymer translocation through 
membranes. 
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IX. CONCLUSIONS 

In this paper, the statistical mechanics of a long flexi- 
ble polymer dragged into a cylindrical tube with repulsive 
wall under good solvent conditions is studied, both via 
a scaling theory and by Monte Carlo simulations, using 
the PERM algorithm that allows the successful study of 
very long chains. It is shown that in the limit of infinite 
chain length an entropically driven abrupt transition oc- 
curs when the distance of the chain end from the tube en- 
trance, x, is used as a control parameter. Experimentally, 
this situation could be realized e.g. when a nanoparti- 
cle is attached to this chain end and the position of this 
particle is controlled externally by a laser tweezer. This 
has to be distinguished from the case when the chain is 
driven by applying a constant force (e.g. electrophoret- 
ically). It is shown that a critical value x c exists, such 
that for x < x c a finite fraction of monomers is outside 
of the tube in a mushroom-like configuration, and the 
remaining Aj mp < N monomers form the " stem" , a one- 
dimensional stretched string of blobs, while for x > x c the 
"crown" of this flower- like conformation of the polymer 
has disappeared, and all N monomers have been sucked 
into the tube to become part of the " stem" , the string of 
blobs. This transition for N — > oo is strongly discontin- 
uous, since at the transition Ni mp /N jumps from about 
3/4 to unity in our model. We construct a suitable or- 
der parameter for this transition and use scaling ideas 
to formulate the Landau free energy branches of both 
states that compete near the transition with each other. 
The Monte Carlo results confirm the general picture of 
these two phases and allow to estimate the undetermined 
prefactors of the Landau theory description. The Monte 
Carlo results also allow to quantify the extent of finite 
size rounding of the transition that inevitably occurs as 
a consequence of the finiteness of the chain length. Ar- 
guments are presented that for cases of physical interest 
(such as DNA in artificial nanopores) it is rather likely 
that this transition is affected by hysteresis, and esti- 
mates for the lifetime of metastable states are given. 

Interestingly, no transition is predicted for a reverse 
process of chains release if the constraint fixing the chain 
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Appendix: A dynamic picture of polymer escape 
from a tube 

We have demonstrated that descriptions of the process 
of polymer chain de-confinement upon release may differ 
depending on the statistical ensemble. The choice of an 
appropriate ensemble depends, in turn, on which of the 
variable: the distant end position, x, or the number of 
imprisoned segments, N- lmp , is the slowest. In a more 
consistent approach, both variables are treated on equal 
basis. We define the Landau free energy as a function 
of two independent variables, &(x, -/Vj mp ) (assuming, as 
usual, that all other internal degrees of freedom equili- 
brate much quicker) . We have all the necessary informa- 
tion to define $(x, Ni mp ) at hand. For a fully confined 
chain with A; mp = N and far enough inside the tube, 

x > xip , the chain is fully relaxed and its free energy 
is given by Fi mp , see eq 1131 For a fully confined chain 

with x < Xsp , variable x has the meaning of the end- 
to-end distance and the expression derived for the flower 
branch of the Landau free energy, eqUU applies. The 
same expression applies for any x provided the chain is 
only partially confined. This eventually gives 

*/ Al \ { 5-40, forJVimp = N, x > 0.9ADn b 
nx ' JVimpJ = I 1.48-^ + 0.67m 5 - 1 + 1.98m" 1 ) , 

Ub [ otherwise 

(25) 

Here rib = N(D/a) x l v is the number of blobs in a fully 
confined chain, and u is the ratio of two reduced quanti- 
ties u = {x/(Dn b ))/(N imp /N). 
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and 



d<S>(x,N ilnp ) 



, respectively, and the coupled 
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<x>/(n b D) 



FIG. 14: Contour plot of Landau free energy landscape as a 
function of two independent variables $(i,JVi mp ). The cor- 
ridor of the fully confined relaxed chain in the upper right 
corner is artificially broadened for a clearer visual picture. 
Also shown are the parametric trajectories < Ni mp (t) > vs. 
< x(t) > of the de-confinement process for two limiting cases 
D x « D n and D x » D n (solid lines) and for isotropic 
diffusion D x — D n (dashed line). 



A contour plot showing the landscape of &(x, JVi mp ) is 
presented in Figure Q31 States corresponding to a fully 
confined relaxed chain are depicted by a narrow corridor 
in the upper right part of the landscape. The lower left 
corner represents a completely free chain with the lowest 
possible free energy. The corridor is generally separated 
from the sloping landscape by a barrier, except for an 

opening in the vicinity of x — xip^ 1 . 

Full dynamic evolution in the (a;,A^ mp ) configuration 
space will be governed by a Fokker-Planck equation with 
Landau free energy playing the role of effective poten- 
tial. However, if one avoids the detailed description of 
barrier-crossing events, the problem is simplified consid- 
erably. The dynamics of averaged quantities < x(t) > 
and < Ni mp (t) > is then driven by thermodynamic forces 



dx a,llu dN in 
non-linear equations of motion are given as 



and 



d < x(t) > 
dt 



d < N imp (t) > 
dt 



-D 



dx 



■D r . 



9$ 



dN h 



(26) 



(27) 



where D x and D n are the effective diffusion coefficients 
along the x and N- lmp coordinates, correspondingly. Ge- 
ometrically, this evolution is a diffusive slide along the 
slopes of the Landau free energy in a two-dimensional 
configuration space. A detailed discussion of the diffusion 
coefficients and the physics behind them is well beyond 
the scope of this paper. 

Here we briefly present the result of a dynamic anal- 
ysis for three different scenarios. First we analyze the 
two limiting cases discussed in Section 8. Equations of 
motion were solved numerically for very slow x dynam- 
ics with D x = lO~ 3 D n and the curve of < Ni mp (t) > 
vs. < x(t) > parametrically defined by t is shown in Fig- 
ure[T31 It is clear that this coincides with the quasi-static 
trajectory predicted in the a;-ensemble. Another curve 
was obtained from a numerical solution assuming slow 
-^Vimp dynamics with D n — 1Q~^D X . This coincides pre- 
cisely with the TVjmp-ensemble description as evidenced 
by comparing Figures IT51 and [hH Finally, the parametric 
curve for the case of isotropic diffusion D x = D n is also 
presented. This is characterized by a relatively rapid 
initial growth of the ejected part of the chain, the dis- 
tant end position evolution catching up with some delay. 
Eventually, the trajectory slide down the valley along its 
geometrical bottom line. The shapes of all three tra- 
jectories are insensitive to assumptions made about the 
change in diffusion coefficients during the dc-confincmcnt 
process as long as their ratio is kept fixed. 
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